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Quadrotor by nature is a very unstable system and flying it without any 
feedback control algorithm is deemed impossible. However, before designing 
the control system, system identification need to be conducted as the 
accuracy of the control system depends highly on the accuracy of the model. 
Therefore, this paper explained the design of the quadrotor model with an 
“X” configuration using the Euler-Newton model. Two types of test rig were 
designed to measure the thrust coefficient, torque coefficient and throttle 
command relation parameter needed in the model. Other parameter such as 
moment of inertia was also being measured by separating the quad rotor 
model into several sections: Motors, Electronics Speed Controllers (ESC) 
and Central Hub. All parameters needed in the designed quad rotor model 
has been successfully identified by measuring the parameters using the 
custom-built quad rotor and test rigs. The parameters found in this paper will 
be used in designing the control system for the quadrotor. 
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1. INTRODUCTION 

Recently, research community tend to focus on unmanned aerial vehicles (UAV) as it has a great 
potential and has been applied in many applications such as security surveillance, disaster management, crop 
monitoring, on-board imaging and educational purposes [l]-[4]. Quadrotors have become very popular due to 
their small size, simplicity of mechanics and good maneuverability [5]. Even though quadrotors can bring a 
whole lot of new application cluster, the research has been stalled until recently. The quadrotor is 
significantly under actuated since it only has four rotors (actuator) to completely control 6 degrees of 
freedoms. The degrees of freedom to be controlled are three rotational, Yaw (0), Pitch (6 ) and Roll (0) and 

3 translations, X, Y, Z. 

The modeling and control system development of the quadrotor has been covered by many scholars 
in their work with varying complexity [6]-[8]. For instance, in [9], the author modeled the quadrotor using 
Newton-Euler model and the modeling includes the fundamental dynamic model of the motors, while [10] 
have modeled the quadrotor only for z-axis maneuver empirically and of first order only. The author of [11] 
also modeled the quadrotor using Newton-Euler model, however modeled the rotor dynamics empirically 
which makes the model hardware specific. Other scholars used Euler-Newton model to model their quadrotor 
dynamics also include [12], [13] and [14]. The Euler-Newton method is fundamentally accurate and is 
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suitable to be used in simulation analysis. When it comes to actual implementation, the process becomes over 
too difficult practically but, some parameter identification can be simplified by estimating these parameters 
empirically as discussed by [15]. 

It is almost impossible to define the quadrotor parameters based on passport data components that 
are usually included in its composition, as the detailed information about their product rarely provided by 
manufacturers. In addition, the dynamic characteristics of the quadrotor system depends on the union 
characteristics of all its elements especially to those system that are custom built. Therefore, the special 
measuring test rigs [16]-[20] are used to solve this problem by providing the measurement of thrust force and 
the aerodynamic torque that are hard to be measured. It focused mainly on the identification parameters of 
fixed step propeller drives system. Therefore, test rig is a system that are required to allow quadrotor 
parameter measurement, as a physical plan and it is being held without linear displacement and usually have 
three degrees of freedom. 

In this paper, a parameter identification process was conducted by measuring the value of mass 
moment of inertia (J), thrust coefficient (6V), torque coefficient (C Q ), throttle command relation (C R ) and 
motor time constant. Two custom test rigs were built to measure the C T , C Q and C R parameter based on the 
principle of thrust and torque coefficient. The parameter found in this paper will be used to represent a 
quadrotor system with an “X” configuration. The parameter identification process was conducted by taking 
the real quadrotor parameter measurement in the lab. The quadrotor model used in this paper will be 
explained in Section 2 and parameter identification process will be described in detail in Section 3. 

2. QUADROTOR MODEL 

In this paper, the configuration chosen for the quad rotor is the “X” configuration. By using the 
vector coordinate notation, the translational position and angular position of the quad rotor with respect to the 
Earth frame can be defined in Equation 1 and 2 respectively. 


11 = [x y z] T 

(1) 

5 = [4>erp? 

(2) 


The acceleration of the center of mass of the rigid body quad rotor is based on the forces and accelerations 
acting on the body are determined by the velocity state equation as described in Equation 3. Where, v is 
linear acceleration of the center of mass in the body frame with respect to the earth frame, m is the total mass 
of the quad rotor, g is the gravity, C is the Euler Angles rotation matrix required to translate the gravity to act 
in the body frame, H is the cross-product matrix for rotational velocity required to map rotational velocity to 
the time derivative of the angles. 




U' 

V 

■W- 


( 3 ) 


When considering the motion of the quad rotor, another state equation that need to be considered is the 
angular velocity state equation as in Equation 4. This equation is based on Eulef s equation for rigid body 
dynamics. This equation describes the change in roll (P), pitch (Q) and yaw(R) rates caused by the inertia, 
angular velocity, and the moments applied by the motor-propeller system. Where H is the cross-product 
matrix for rotational velocity and (x) b is the rotational velocity of the quad rotor body within the body frame. 


cb = (jy 1 [M- ajco b ] = 


P 

Q 
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( 4 ) 


Thrust and aerodynamics applies force on the quad rotor on the positive z-axis as given by Equation 5. 


F = 


0 

0 

c T (gti + gt| + gt| + ml). 


( 5 ) 


The relationship matrix between the thrust and torque with motor RPM can be constructed in Equation 6, 
where d is the length of the arm. 
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( 6 ) 


To determine the parameters, an actual hardware is needed to be built for which the parameters of the model 
are to be determined empirically. In terms of Microcontroller unit, Arduino DUE was used. For gyroscope 
and accelerometer unit, an L3G4200D and ADXL335 from Sparkfun was used. To drive the EMAX 1047 
propeller, EMAX XA2212 motor was used. For ESC, Hobbysky model was used. In terms of frame, DJI 
F450 and Hobbyking HK-T4A was used as a remote controller and EMAX battery was used. 


3. PARAMETER IDENTIFICATION 

In this section, the equations related to the dynamic model of the quad rotor were discussed and of 
which the parameters to be identified from actual physical components are ] h , m, c T , c Q) c R , b and tc. These 
parameters need to be determined empirically and will be discussed in detail throughout this section. 

3.1. Mass Moment of Inertia, J 

To determine the moment of inertia, the quad rotor vehicle was separated into different components. 
The model of each components was represented by a simplified geometric shape of constant internal density. 
Then, the weight of each component was measured. The moment of inertia of each component about the x, y, 
and z axes of the vehicle is determined by using the parallel-axis theorem. The total moment of inertia matrix 
of the vehicle, is the sum the inertias for every component about each axis. The component that needs to be 
modelled are the motor, the ESC, arms and the central hub. The parallel axis theorem is defined in Equation 
7. 


J parallel-axis = Jcom + mr 2 (7) 

Where, ] C0M is the moment of inertia of an individual component along its own axis parallel to the axis it is 
to be moved, m is the mass of individual component and r is the perpendicular distance between the parallel 
axes. The modeling of each individual component of quad rotor and their mass moment of inertia along 
corresponding axis is calculated in the next section. 

3.1.1. Motors: Solid Cylinders 

The motors were modelled as a solid cylinder to determine the mass moment of inertia for all four 
motors. Figure 1 shows the quadrotor structure that indicates the dimensions that need to be considered to 
estimate the mass moment of inertia for motors. Equation 8 shows the equation of the mass moment of inertia 
for x and y axes and Equation 9 for z axes respectively. 


m - mass of one motor 
d m — distance from motor to C.O.M. 
h — height of motor above arms 
r— radius of motor 



Figure 1. Dimensions to be considered for mass moment of inertia calculation of the motors 


Jx - Jy - 2 [^mr 2 + ^rn/i 2 ] + 2 [^rnr 2 +^mh 2 + md m 2 ] (8) 

Jz = 4[imr 2 +md m 2 ] (9) 
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Based on the specifications that are used for the motors in this experiment, the mass (ra), dimensions (d m ), 
height (/i) and radius (r) was found as m = 57 g, d m - 26.035 cm, h= 3.0 cm and r= 1.35 cm. Using the 
dimensions measured, the moment of inertia of motor corresponding to each axis is ] x = ] y = 7.8059 X 
10" 3 and ] z = 0.015475. 

3.1.2. ESC’s Thin Flat Plates 

The ESC’s were modeled as thin flat plates and the mass moment of inertia for all four motors is 
found based on Figure 2, Equation 10 and 11. Based on the ESCs used in this project, the mass and the 
dimensions of the ESCs were measured to be as m = 22 g, a = 4.5 cm, b = 2.6 cm, d s = 11.05 cm. Using the 
dimensions measured, the moment of inertia of ESCs corresponding to each axis is ] x ~]y — 5.47154 X 
10" 4 and J Z} = 2.17778 x 10" 5 . 


Motor 4 


Z 


m — mass of one ESC 

a - width of ESC 

b — length of ESC 

d s — distance from ESC to C.O.M. 



Figure 2. Dimensions to be considered for mass moment of inertia calculation of the ESCs. 


] X =jy = 2 [^ma 2 ] + 2 mb 2 + md s 2 ] 

(10) 

Jz = 4 m(a 2 + b 2 ) + md s 2 ] 

(ID 


3.1.3. Central HUB: Solid Cylinder 

The central hub was modeled as a solid cylinder and the mass moment of inertia for the central hub 
was found based on Figure 3, Equations 12 and 13. 


m — mass of central hub 
r— radius of central hub 



Figure 3. Dimensions to be considered for mass moment of inertia calculation of the Central HUB. 


Jx=Jy=[lmr 2 +±mH 2 ] 

(12) 

Jz = [lmr 2 ] 

(13) 
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Based on the quad rotor, the mass and the dimensions of the Central HUB were measured to be as m = 335 g, 
r = 6.2 cm, H = 8.0 cm. Using the dimensions measured, the moment of inertia of the Central HUB 
corresponding to each axis is as J x = J y = 5.006 x 10 -4 and J z = 6.4387 x 10 -4 . 

3.1.4. Arms: Long Cylindrical Rods 

To measure the moment of inertia for the arms, a cylindrical rods configuration was selected. The 
mass moment of inertia for the arms were found based on Figure 4, Equations 14 and 15. 


m — mass of one arm 
r — radius of arm 



Figure 4. Dimensions to be considered for mass moment of inertia calculation of the Arms 


J x =J y = 2 mr 2 ] + 2 [^mr 2 + j mL 2 + md/] (14) 

J z = 4 [^mr 2 + ±mL 2 + md/\ (15) 


Based on the Frame, the mass and the dimensions of the Arms were measured to be as m = 57 g, r = 0.5 cm, 
L = 22 cm, d a - 3.5 cm. Using the dimensions measured, the moment of inertia of the Arms corresponding to 
each axis is ] x = J y = 1.98 X 10 -3 and J z = 3.959 x 10 -3 . Finally, by summing ] x of all individual 
components, J y of all individual components and the ] z of all individual components, the moment of inertia 
matrix is populated in Equation 16. 


7 = 


0.0108 0 0 

0 0.0108 0 

. 0 0 0 . 0212 . 


(16) 


3.2. The Thrust Coefficient, C T 

The thrust coefficient can be found by determining the gradient of the thrust versus (RPM) 2 graph. 
The principle of measurement of thrust can be seen as in Figure 5. The mechanism is pivoted at point A and 
as the thrust is generated by the motor, it can be measured at point B. 


Thrust of Motor 



Figure 5. Principle of thrust measurement 


A test rig was built to visualize the principle of measurement in Figure 5. To measure the thrust 
coefficient, a test rig shown in Figure 6(a) was built. To capture the RPM value, an electronics system was 
designed as shown in Figure 6(b). The RPM value was simply captured by using phototransistor, which 
detects the blade crossing, and RPM value was determined by measuring the time between the consecutive 
blade crossings. This data was then fed into an Arduino based microcontroller and then imported to 
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Microsoft Excel to be plotted into a graph of thrust versus (RPM) 2 . To capture the thrust, a simple household 
digital weight was used. 



Figure 6. (a) Test rig to measure Thrust Coefficient, (b) Electronics system to measure RPM 


To find the thrust coefficient, the Thrust value vs RPM 2 graph was drawn and the best-fit gradient was then 
determined. Figure 7 shows the data captured on the graph of thrust vs (RPM) 2 . Then, a linear approximation 
was made and finally the thrust coefficient was found to be 2' 7 . 
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Figure 7. Graph of Newton vs (RPM) 2 to determine the thrust coefficient. 


3.3. The Torque Coefficient, C Q 

The principle of measurement of Torque Coefficient was shown in Figure 8. 



PI A PT B 


Figure 8. Principle of torque measurement 


To visualize the principle of measurement in Figure 8, another test rig was built as shown in Figure 9. The 
torque coefficient was determined by drawing the best-fit gradient in the graph of Nm vs RPM 2 . 
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Figure 9. Test rig to measure the torque 
coefficient 


Figure 10. Graph of Newton*m vs (RPM)2 to determine 
the torque coefficient 


Figure 10 shows the graph of torque generated by the corresponding RPM 2 . After performing the linear 
approximation, the torque coefficient was found to be 2‘ 9 . 

3.4 Throttle Command Relation, C R 

To measure the throttle command relation, either test rig used to measure C Q or C T parameter can be 
used. To determine the parameter C R and b , the best-fit gradient line was drawn on the graph of RPM vs 
Throttle. 


% Throttle vsRPM 
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STEP RESPONSE OF 20% THROTTLE TO 3034 THROTTLE 



Figure 11. Graph of throttle percentage vs RPM Figure 12. 10% step response to determine the motor 

time constant 


By plotting the graph of throttle percentage versus the RPM and performing the linear approximation as 
shown in Figure 11, C R was found to be 107.65 and b, which is the y intercept of the graph is found to be 
1302.7. 

3.5. Step Input to Determine the Motor Time Constant 

The time constant of the motor is used to model the dynamics of the motor in simulation. To 
determine the motor constant, a step of 10% throttle is applied and the response was captured. Then, the time 
constant was determined at T63% from the step response graph. Figure 12 shows the step response graph and 
the constant was found to be 0.07 seconds as TO = 25.385 s and T63% = 25.455 s. 


4. CONCLUSION 

In this paper, an approach was taken to determine the parameters required for populating the 
mathematical model. Two custom test rigs were built by visualizing the principle of thrust and torque to 
measure the thrust and torque coefficient. The parameters were determined as ] x / y = 0.0108, ] z = 0.0212, 
C T = 2 -7 , C Q = 2~ 9 , C R = 107.65 , b = 1302.7 and t c = 0.075. All parameters that were found in this 
paper will be substituted in the mathematical model for controller design in future works. The controller 
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based on the parameters that were found in this paper will be designed and the simulation and experimental 
test will be expected to be conducted to evaluate the performance of the designed controller. 
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